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Abstract 



Q 
U 

I We present exact analytical results for the statistics of nonlinear coupled oscillators 

under the influence of additive white noise. We suggest a perturbative approach 
for analysing the statistics of such systems under the action of a determanistic 
perturbation, based on the exact expressions for probability density functions for 
^ ' noise-driven oscillators. Using our perturbation technique we show that our results 

00 . can be applied to studying the optical signal propagation in noisy fibres at (nearly) 

zero dispersion as well as to weakly nonlinear lattice models with additive noise. 
' The approach proposed can account for a wide spectrum of physically meaningful 

O ■ perturbations and is applicable to the case of large noise strength. 

o 

■ Key words: Stochastic dynamics, nonlinear oscillators, Fokker-Planck equation, 
nonlinear optics 

■ PACS: 05.10.Gg, 42.81.Dp 
> ' 

X 



1 Introduction 



Weakly nonlinear coupled systems (discrete self-trapping model, discrete non- 
linear Shrodinger equation, etc.), belong to the universal, widely applicable, 
highly illustrative and thoroughly studied models of nonlinear physics. These 
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models have applications to molecular crystals, molecular dynamics, nonlin- 
ear optics, biomolecular dynamics, and so on (see e.g. [1,2,3,4] and references 
therein). It was shown that they exhibit unusual dynamical phenomena per- 
taining only to the nonhnear discrete systems, hke existence of intrinsic lo- 
calised modes etc. [5]. Apart from the dynamical behaviour the influence of 
random force and statistical properties of nonlinear lattices are also subjects 
of keen interest [6,7,8,9]. The complex dynamics of such systems is determined 
by an interplay between such factors as randomness, discreetness and nonlin- 
earity. The results concerning the statistics of coupled nonhnear oscillators 
driven by Gaussian white noise, given below, can be straightforwardly applied 
to the class of problems which are naturally discrete, e.g. for the noisy self 
trapping model and its modifications. The second (but not less important) 
purpose of the current paper is to show that the stochastic discrete models 
with flnite number of degrees of freedom can also be applied to the regularised 
stochastic continuous systems and can have a much broader fleld of applica- 
tion. More specifically, in the spirit of the original idea of Mecozzi [10,11], we 
show that the model of randomly driven coupled nonlinear oscillators can be 
used for studying the statistics of signal propagation in a nonlinear optical 
fibre with inline noise sources. 

A universal model describing the propagation of the wave envelope in a weakly 
nonlinear dispersive media (which include optical fibres) is the Nonlinear 
Shrodinger Equation (NLSE) and its scalar and vector modifications [12,13,14]. 
The need for the profound exploration of noise-stimulated fiuctuations of sig- 
nal parameters in an optical fibre stems from the great practical importance 
of such study for immediate technical purposes [12,13]. The pulse propaga- 
tion inside an optical fibre with inline amplifiers is described by the perturbed 
NLSE with additive noise, which accounts for amplifier spontaneous emission 
(ASE)[12]. It provides an illustrative example of a complicated continuous 
noisy system where the exact results concerning the statistics of the signal 
are still very scarce owing to the complexity of the problem (see for instance, 
review article [15] and also [16,17] for recent results). However if we assume 
that the noise is delta correlated in both time and space a natural regularisa- 
tion occurs, which transforms our continuous system into a system of coupled 
nonlinear oscillators subjected to an external white noise. The statistical prop- 
erties of such regularised discrete system are much easier to derive than those 
of the original continuous systems. In particular one can use the Fokker-Planck 
equation (FPE) approach to obtain the evolution of the probability density 
function (PDF) for all discrete components of the signal. But even the dis- 
crete system corresponding to noisy NLSE and its modifications is still very 
complicated because of the coupling between the oscillators, which is due to 
the dispersion. Therefore Mecozzi [10,11] considered propagation of an optical 
pulse affected by joint action of Kerr nonlinearity and additive white Gaussian 
noise (WGN) at zero dispersion, i.e. the dispersive (second time derivative) 
term in the NLSE was dropped. The resulting effective dispersion free equa- 
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tion (after the regularisation procedure) formally coincides with the dynamical 
equation for a single nonlinear oscillator driven by a noisy external force. 

In the present paper wc treat a zero dispersion stochastic system (a single 
noise-driven nonlinear oscillator) and its straightforward generalisations (non- 
linearly coupled noisy oscillators) as basic systems and consider various small 
deterministic perturbations which can be of quite arbitrary nature: higher or- 
der nonlincarities, small residual dispersion, etc. In particular we can add a 
perturbation which will be the direct discrete analogue of that of the corre- 
sponding continuous optical system. Other continuous models yielding nonlin- 
ear discrete equations addressed in the current paper, include higher dimen- 
sional systems of nonhnearly coupled NLSEs. Two nonlinearly coupled oscil- 
lators with particular coupling coefficients correspond to well known Manakov 
equations (see e.g. [18]) with zero dispersion driven by white noise. This con- 
tinuous model describes the nonlinear pulse propagation in a noisy birefringent 
fibre at nearly zero average dispersion. Systems of more than three coupled 
NLSEs are also of physical relevance [19]. In addition to optical communica- 
tions, in the context of biophysics the case of three nonlinearly coupled NLSEs 
can be applied for studying the propagation of solitons along the three spines 
of an alpha helix in protein [20]. Systems involving higher numbers of coupled 
NLSEs have applications in the theory of optical soliton wavelength division 
multiplexing [21], multichannel biparallel- wavelength optical fibre networks 
[22] , and so on. Hence we believe that the study of the stochastic properties of 
such models (which, to the best of our knowledge, is a relatively undeveloped 
field) is also of great interest and may have a significant practical outcome. 

As mentioned earlier, we describe the system statistics in terms of the proba- 
bility density function for the discrete variables and the evolution of the PDF 
is governed by the FPE [23] . In the case of an unperturbed system of noise- 
driven nonlinear oscillators (which is a discrete analogue of a system of zero 
dispersion noisy NLSEs), the corresponding FPE can be solved analytically. 
The main idea of our approach is to apply the perturbation theory directly to 
the FPE rather than to the initial stochastic system, using analytical results 
for the base system. The advantage of such a method is that it allows one 
to obtain the PDF for the signal output directly as a series in powers of a 
small parameter e, which is the effective "strength" of the perturbation. The 
perturbation theory is based on the propagator of the unperturbed FPE. Once 
derived, the propagator can be applied to various discrete systems, e.g. regu- 
larised systems of weakly coupled NLSE at zero dispersion, where the couphng 
is considered as a perturbation. This means we can incorporate the effects of 
coupling and consider the statistics of nonlinear lattice models in the so-called 
anticontinuum limit (weak intersite linear coupling). In the context of optical 
applications we are able to consider real polarisation mode dispersion (PMD) 
systems as well as systems with wavelength division-multiplexing (WDM) [13]. 



3 



The paper is organised in the foUowing way. In section 2 we introduce basic 
models of one and more noise-driven coupled nonlinear oscillators and write 
down the Fokker-Planck equation for each model. After that we consider noisy 
scalar and vector NLSE, describing the propagation of a signal in the optical 
fibre, and show that after the regularisation (i.e. after the introduction of a 
discrete time variable) the system becomes equivalent to one of the "base" per- 
turbed systems. We also mention the stochastic discrete self-trapping model 
as an example of naturally discrete system which can be analysed perturba- 
tively. In section 3 we proceed to find the propagators of FPEs derived for 
noisy oscillator systems. Using the obtained propagators in section 4 we ad- 
vance to build a perturbative expansion pf the PDF for different types of 
perturbations. In this section we consider non-dispersive perturbation of one 
and more coupled zero dispersion NLSE and provide explicit expressions for 
several typical example systems. In section 5 we move to more complicated 
systems, i.e. to the scalar or vector NLSEs, where second dispersion is treated 
as a perturbation. In the conclusion we summarise the results and outline the 
key features and perspectives for the approach proposed. 



2 Basic models, equations and regularisation procedure 

In first two subsections we write down the explicit form for the FPEs attached 
to the system of stochastic equations, governing the evolution of oscillator 
statistics. Then we explain how these systems can be employed in the context 
of optical propagation and nonlinear lattices. 

2.1 Single noisy nonlinear oscillator 

By adding a white noise term to the (dimensionless) dynamical equation for 
a nonlinear Kerr oscillator one gains the Langevin equation for the complex 
field u = X + iy: 



In the mechanical interpretation x is a position and y is a velocity of the 
nonlinear oscillator or x and y are the components of a torque for the nonlinear 
rotator; has a meaning of time and ri{Q = fJiiO + ^V^iC) is the complex white 
Gaussian noise with the following correlation properties: 



with D being the noise strength. Equivalently one can decompose the real and 
imaginary parts of (1) and obtain a system of two coupled Langevin equations: 




i\u\\ + r){Q . 




<Vi{C)r]j{C') >= 2DS,jS{C-C), 



(2) 
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^ = + 77,(0, - {x' + y') x + r„{0 . (3) 

Introducing the 2D vectors q — {x, y} and rj — {rji, 772} we can write equation 
(3) in general vector form: 

^^fiq)+viO- (4) 

Here / is a deterministic "advection vector" (note that sometimes in the 
hterature it is written with a minus sign). This is the canonical form of the 
Langevin equation with additive WGN. We will be interested in the statistics 
of the solution of Eq.(4) for a given (. Namely, we will seek the conditional 
PDF of the vector q, P{q, (), provided that at the initial moment (( = 0) the 
PDF is given. This function (see [23]) obeys the following second order partial 
differential equation: 

^ = i:{-di[Mq)P]+DdfP}, (5) 

which is a FPE attached to system (4) with the additive white noise, given 
by (2). Also one has to impose the normalisation condition JqPdq = 1. A 
propagator is a special solution of (5) with the initial condition P{q, 0) = 
S{q — q°), i.e. it corresponds to the conditional PDF provided that at the 
initial moment the system is in a deterministic state q^. The explicit form of 
the FPE for single oscillator system (1) is: 

(6) 

The boundary conditions for (6) are chosen as follows: the PDF, P{q, (), must 
have no singularities and should decrease rapidly as |q| — > 00 to provide the 
normalisation. 



2.2 System of nonlinearly coupled oscillators 

Next we study the following system of M coupled nonlinear oscillators: 




J At 

1 • X — ^ I 1 2 1 / >\ 

^iui2_^aii\ui\ +7] (C) 
•^s i=i 



dUM . I |2 I M/^\ 

-^^iUM2^aMi\ui\ +ri {Q 

i=i 



(7) 
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Here again Ui = Xi+iyi, arc the complex fields, aij arc real constant coefficients 
and the correlation properties for complex WGNs rj'' = 7]\ + ir^g arc (here and 
further on the subindexes of complex WGNs, "1" and "2" , imply the real and 
imaginary parts correspondingly): 



(8) 



(here no summation over the repeated index k implied) . To be more general we 
suppose that WGNs r]^{C,) may have (possibly) different intensities Dk. Again 
we can introduce 2M-dimcnsional real vectors q = {xi, . . . , xm, Vi, ■ ■ ■ Vm} and 
T] = {rjl, . . . 7]^ , vh ■ ■ ■ v¥}i after singling out real and imaginary parts in 
(7) arrive at the generic Langevin equation (4). The FPE for this system has 
the form 



dP 



Ml M 



i=l 



d_ 



d_ 

'dxi 



(9) 



Here Aj = We shall call the number M of nonhnearly coupled 

equations in (7) the dimensionality of the system (the number of independent 
variables in the FPE (9) is then 2M). 



2.3 Regularisation procedure for optical pulse propagation: discrete time 



The propagation of a complex light envelope u{(,t) in a noisy optical fibre line 
with so-called weak dispersion management (see Refs. [12,24]) is described by 
the stochastic NLSE: 

(Note that here C plays the role of coordinate in contrast to the mechanical 
interpretation of oscillators and "virtual" time is t .) In equation (10) the 
stochastic term ri{t, C) represents a complex WGN. The following correlation 
properties are assumed: 

<77(i,C)>=<77(i,CM^',C')>=0, 

<r?*(t,CMt',O>=2G'5(C-C')<^(^-0- ^ ^ 

Here G is a normalised noise strength. Equation (10) describes the path- 
averaged model of the dispersion managed optical fibre communication system 
and, it was shown (see [24] and references therein for details) that this model 
is applicable when the so-called strength of the map, proportional to the lo- 
cal fibre dispersion, to the dispersion oscillation period and to the squared 
reciprocal pulse width, is small. This strength parameter characterises the ef- 
fect of the variation of local dispersion. Parameter /9 entering Eq.(lO) is just 
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a path-averaged dispersion coefficient. A zero dispersion NLSE occurs when 
the path-averaged dispersion coefficient is equal to zero, [3 = 0. The statis- 
tics of system (10) at zero dispersion point were first studied in [10,11] (see 
also [25]). In these Refs. an analytical formula for the conditional PDF of the 
output signal u{t, L) given the deterministic input signal uo{t, 0) was obtained. 

Here we would like to consider the zero-dispersion NLSE assuming also that 
there is an addition corresponding to a deterministic perturbation: 

— = i\u\\ + ri{t,C)+eN[u,...]. (12) 

Here £ <S 1 is intended and let N[...] be any sort of perturbation of the initial 

zero-dispersion NLSE: it may be either a "nonlocal" operator involving the 
time d/dt derivatives N[u, u'^, u'l^, . . .], or some kind of a "local" operator, N[u]. 
Note that complex WGN defined by expression (11) is equivalent to the two 
independent real WGN, rj{t, () = rii{t, () + iri2{t, (), with 

< 771 (t, C) >= < V2{t, > = < Vl{t, ()V2{f': C) >= , 

< vi{t, Oviit', CO >= < V2{t, 0v2{t', CO >= G6{c - C)S{t -t'). ^ ^ 

Because the white noise is delta correlated in both time and space, it has an 
infinite average power < |r7(t, C)P >• This is the consequence of the idealised 
character of the white noise. In practice, real noise always has small but finite 
correlation length Tc and correlation time r^- For ASE in optical fibre trans- 
mission links the correlation length is in the order of the average amplifier 
spacing L^, and correlation time is inversely proportional to the bandwidth 
of the noise: Tc ~ B^^ . In what follows we will consider our noise still delta 
correlated in space. However to make the correspondence between discrete and 
continuous systems exphcit and also to get physically meaningful results we 
move to the context of discrete signals and introduce discrete time variable t. 
We will consider the case where both signal and noise are entirely contained to 
the time interval of length T. Then, to sample the signal effectively, we need a 
finite number of discrete samples, separated by an interval At ~ Tc ^ T. If we 
put At = -B~^, the number of samples is iV = BT. The signal is represented 
as A^-dimensional complex vector, u = a? -|- it/, in the sample space: 

«(t,C)^w(C), «^(C)=«(^A^,C), i = l,...N. (14) 

Analogously one can work in the real 2A'^-dimensional space of real and imag- 
inary parts of vector u: q = {xi, . . . , xjy, yi, . . . , i/n}- After the regularisation 
Eq.(12) is equivalent to the system of equations: 

-^ = i\ui\\i + ri\C) + eNi[u], (15) 

Vector t;(C) = ''7i(C) + i^2(C) is a regularised white noise with the same prop- 
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erties as given by (13) but with regularised temporal delta functions: 

< vmvm >=< vmviiC) >= 205,, 6{c - a, m 

where we set D = BG. The operator TVj is discretised in time as well. 

We call the equation, governing the unperturbed stochastic evolution of each 
separate sampled component, the base system. For all cases considered in this 
paper, a base system will be either of the form Eq.(l) or Eq.(7); the corre- 
sponding FPE for each base system has either a form (6) or a more general 
look (9), depending on the dimensionality of the base system. 

One can see that each sampling component in Eq.(15) obeys the equation 
governing the noisy dynamics of a single nonlinear oscillator, Eq.(l): hence 
the dimensionality of the corresponding base system is one. The complete set 
of equations for all components of the sampled signal has a form of perturbed 
Eq.(7) with M — N and an identity matrix of a-coefficients: = 5ij. As was 
noted before, we divide the perturbations in equation (15) into two categories: 
"local" , which depend on the value of the signal at the current moment of time 
only, i.e. Ni[a\ = Ni[ui], and "nonlocal" (as, for instance, dispersion, involving 
second or higher time derivatives of the signal) , which comprise the field values 
at other sampling points. 

The same regularisation procedure can be performed with other models rele- 
vant to the noisy signal propagation. For example we can consider the propaga- 
tion in the birefringent fibres where instead of NLSE, under certain conditions, 
we have a set of Manakov equations [18] perturbed by noise: 

dv B&^U . /, ,9 , ,9\ 9/^ X 

9C^"'^9^ + ^H'''' )+'^ 

Here u and v stand for complex components of the left and right polarised 
waves and rf{C,, t) are the complex WGNs. Without noise addition this system 
describes the processes of the two-dimensional stationary self-focusing and 
one-dimensional auto-modulation of electromagnetic waves with arbitrary po- 
larisation. Neglecting the dispersion and adding perturbation, after the regu- 
larisation one obviously gains: 

i(|^/i|' + \vi\'')ui + riliC) + eNl[u, v] , 

(18) 

i{\ui\^ + 1^;^!^)^;^ + riiO + eNl[u, v] . 

In this case, for each separate sampled pair [wj,!'!] (if £ is again set to zero), 
we arrive at the system of two nonlinear ly coupled oscillators. This is just a 



duj 

dC 
duj 
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special case of Eq.(7), where the number of dimensions of the base system 
is two. The whole system for K = 2N sampled components of the signal 
[N is, as previously, the number of the sampling time-points) has the form of 
Eq.(7) as well. However, as the base system has a dimensionahty two, here the 
coupling matrix aij is no longer an identity matrix but have a (2 x 2)-block 
diagonal form. The "local" perturbations for such systems are those which 
satisfy: Ni[u,v] = Ni[ui,Vi]. 

In the same way one could deal with the systems involving higher number M 
of coupled NLSEs: the number of equations would simply define the dimen- 
sionality M of the corresponding base system. The regularised system for the 
whole set of time-sampled signal components would possess the matrix Ojj in 
a block- diagonal form, where each block would be an (M x M) square matrix. 

2.4 Nonlinear lattice model 

The dynamics of a weakly nonhnear lattice is described by the set of differential- 
difference equations (the discrete self-trapping equation (DSTE), see Refs.[l,2,4]): 

^ = i\uifui + ieY^rriijUj + ri\C) , (19) 

where we have also added the white noise terms '>]''{() at each site (with the 
same statistical properties as in (16)). In Eq.(19) rriij is a coupling matrix and 
index j ranges over the whole ID lattice containing N sites. The generahsation 
for lattices of higher dimensions is straightforward. In physical applications rriij 
is real and symmetric; in the case of tridiagonal matrix, i.e. for the nearest 
neighbor coupling, we recover a particular realisation of DSTE - the discrete 
nonlinear Shrodinger equation [3]. 

Equation (19) provides an example of a naturally discrete system with Kerr 
nonlincarity. If we consider the case of vanishing e, the so-called anticontinuum 
limit, we arrive at the system of uncoupled oscillators, which is obviously the 
unperturbed system (15). Here the dimensionality of the base system is one. 
The perturbation in Eq.(19) provides an example of a nonlocal perturbation, 
since it involves other lattice sites. 



2. 5 The list of the notations 

In this subsection we summarise the notations that will be used through- 
out the paper. The aim of the proposed notations and indexing scheme is to 
demonstrate the general ideology and to underline the similarity of all final ex- 
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pressions: very different physical systems after the regularisation yield similar 
formal expansions for the PDF. 

Tabic 1 



Table of notations and indexing scheme used in the paper 



ISTo^.a i", 1 on 


lVrpaniri0* 


U, V, Ui, Vi 


The sampled complex field variables. The bold symbols 
indicate the whole sets u = {ui, ti2, . . .} etc. The subscript 
indicates the sampling time-point number. 


Xi, Ui and r^, 


Real and imaginary parts and modulus-phase representa- 
tion of the complex field variables. The subscript shows 
the sampling number (in some examples the upper index 
will be used as well). 


(f>, q, q 


Bold symbols denote sets r = {ri,r2,...}, 4> = 
{^^1) ?^>2, • • •}; where indicated explicitly, the subindex 
shows the dimensionality of the set. Vectors q and q are 
the joint sets: q^^ = {xi,X2, ... ,XK,yi,y2, ■■■ ,yK}, and 
qx = {n,r2, . . . ,rK,(l)i,4>2, ■ ■ -Ak}- The superscript "0" 
marks the initial deterministic values. The upper index is 
also used to distinguish between the different sets. 




The coordinate part of the FP-operator. The subscript 
shows the (half )number of independent variables involved. 


yqiSj 


Pm(q; C)) P^{q'i C) are the PDFs for the base system and 

■friT 't'tip wVinlp tSfmrnlpH cjiOTifil Prir rrinvpnipncp wifln fi nib- 

script we shall explicitly mark the dimensionality of the 
underlying base system M, i.e. for Pm we have q = qj^. 
A superscript will be used to mark the whole PDF, and to 
show the total (half) number of the independent variables 
involved. If the dimensionality of the base system is M, and 
the number of sampling points is N , then K = M x N , 
q = qj^. The calligraphic will mark the PDFs of the 
perturbed systems. The tilde above the PDF will mean the 
representation in "polar coordinates", {x,y} {r,(t)}. 
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Table of notations and designations used in the paper (cont'd) 



^.[...],A^[...] 


Ni[. . .] and . .] are the regularised perturbation op- 
erator and its FP-image correspondingly. The explicit 
form of the FP-image is derived using the general rules 
for the construction of the FPE [23] . 




Complex regularised WGN. The superscript is used to 
distmguish between the diirerent independent WGJNs 
entering different equations. 


Qm{- ■ ■)■, Qk{- ■ ■) 


Qm{- • •) is the propagator of the base system FPE with 
the coordinate part Cm- Symbol Qk will mark the 
propagator of the whole FPE of the entire sampled 
system. 



3 Propagator functions of the unperturbed Fokker-Plank equation 

We have just shown that after the regularisation all the models considered in 
the previous section yield same type of system of Langevin equations, which is 
merely a perturbed system of coupled nonlinear oscillators (7). Now consider 
the statistical properties of unperturbed system (7). The FPE attached to 
system (7) is of course Eq.(9). In this section we will calculate a propagator 
for Eq.(9) which will allow us to build an effective perturbative expansion for 
the PDF in the next sections. 

First we note that in some examples (like for instance a scalar regularised 
zero dispersion NLSE (15) with e = 0) the coupling matrix aij is diagonal, 
and the dimensionality of the base system is one. In that case FPE (9) admits 
factorised solutions in the form: 

K 

P'^iQK; C) = ^(^1, ...,XK,yi,... vk; = 11 M^i^ Vi; C) , (20) 

i=l 

where Pi{x, y;() is a, PDF which satisfies unperturbed FPE for a single oscil- 
lator (6). The factorised form of Eq.(20) is a natural manifestation of the fact 
that for diagonal aij all nonlinear oscillators in system (7) become dynamically 
uncoupled and, because of (8), statistically uncorrelated. Therefore in order 
to study the statistics of such a system it is sufficient to study the statistics 
of a single oscillator (1). Bearing this in mind we initially consider a FPE for 
a single oscillator (6) and find its propagator, and then utilise the obtained 
results to find the propagator of generic system (9) by virtue of Eq.(20). 

The same is true for the unperturbed systems for which a base system has a 



11 



higher number of dimensions, say M. In such the couphng matrix aij 

is (M X M)-block-diagonal, and the expression for the total PDF has again a 
factorised form: 

N 

P^(g^; C) = . . . , XK,yi, ■ ■ ■ Vk; C) = U PuiqV. , (21) 

i=l 

where PMiqu'i C) is now a PDF, which satisfies unperturbed FPE for a system 
of nonhnearly coupled oscillators (9) and N is the number of samples so that 
K = N X M. The factorisation in Eq.(21) is possible because each of the N 
block-systems of M nonlinear oscillators is uncoupled from all other remaining 
block-systems. Now to gain the statistics of this system we are to study the 
statistics of coupled oscillators, described by FPE (9). After it is done we 
employ Eq.(21) to find the final answer for the statistics of the entire set of 
samples. 



3.1 Propagator of the Fokker-Plank operator for a single oscillator 



Let us study the equation for a single oscillator (6). Proceeding to the po- 
lar coordinates x = rcos0, y = rsin0 we introduce function Pi(r, = 
Pi(a;(r, 0), ?/(r, 0); C)- Note that function Pi is not the PDF in polar coordi- 
nates since it docs not have the Jacobian r included. For the function Pi we 
now obtain the following equation: 



dPi ^ ~ dPi 
+ CiPi = ^ + 



-DA"" + r' 



d 



Pi =0, 



(22) 



where 




is the Laplace operator in polar coordinates. We define a propagator., Qi, as 
a special solution of (22) (as a function of {C,r, 0}), possessing the following 
properties: 



^i(C,C';0,0';r,r') 



C=C' 



6{(f)-(f)')6{r -r')/r, 



(23) 



gi{C,C';<P,<P';r,r')^o, c<C- 



Note that propagator (23) coincides with the retarded Green function of Eq. 
(22). Since our system is uniform in ( and (p, the propagator depends only 
on the differences A( = ( — (' and A0 = (p — (f)': Qi = Qi{A(; A(f);r,r'). To 
find the propagator we first determine eigenfunctions and eigenvalues of the 
Fokker-Planck operator Ci in (22). The right, ^, and left, ^, eigenfunctions 
are calculated in the Appendix A. Using these eigenfunctions (A. 6), (A. 7), we 
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can write the sought propagator as an eigenvalue expansion: 



n, V 

Using identity (A. 9) the sum over n can be convoluted yielding the following 
result: 

1 °° n/ 



X exp 



2nD sinh(2a,AC) 



2D 



^l'^lUsinh(2a,AC)j- ^^^^ 



Here stands for the z/th order modified Bessel function of the first kind 
and the quantities are defined in the Appendix A. Formula (24) is valid for 
AC > 0. For AC < the propagator, ^i, must be equal to zero due to (23). 
Note that the obtained propagator is symmetric with respect to variables (r, r') 
and is normalised to unity: drr Jq'^ d(j) ^i(AC; — 0'; r, r') = 1. 

Propagator Pi = Qi{( — C'j'P ^ (f)';r,r') (with r{x,y) — \Jx^ + y^, (f){x,y) — 
axctan{y/x) and the same for primed variables) represents the solution of 
equation (6) provided that at the point C — C our variables have deterministic 
values: x = x',y = y' . If we want to apply our results to the regularised 
scalar NLSE with additive white noise, we must use formula (20). Function Pi 
describing the statistics of the unperturbed NLSE at zero dispersion coincides 
with that obtained in [11] with a different approach (see also Ref. [25]). 



3.2 Propagator for the Fokker-Planck operator for M nonlinearly coupled 
oscillators 



In the same way, generalising the derivation of PDF for a single oscillator, we 
can obtain the PDF Puif-, fQ-,4>oi C) for system of M nonlinearly coupled 
oscillators (see Appendix A). FPE (9) in polar coordinates reads as: 



dP, 



M 



+ CmPm — 



dPi 



M 



M 



+ E 



-AA[ 



M 

E 



d 



r^ Uji- 



Pm =0, 



(25) 



where A[ denotes the Laplacian in polar coordinates {r^, (pi}. Similarly to the 
case of one oscillator, we introduce propag ator Gm (AC; (p - (p'; r, r') for Eq. 
(25), which is a solution of (25) for AC > subject to conditions: 



M 



= U^i<l^^-<l>'i)^ir,-r'^/n, 
gM{AC;(j>-4>';r,r')^0, AC<0. 



(26) 
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Again the propagator coincides with the retarded Green function for Eq.(25). 
Using cigcnfunctions (A. 11) and convolution formula (A. 9) we obtain the fol- 
lowing expression: 



e;M(AC;0-c/>';r,r') 



1 



M 



M 



X 



ai 



sinh(2Q;iAC) 



exp 



E E n 

1^1= — OO l'M = — OOi=l 

aiirf + rf) 



(27r) 



2A 



coth(2ajAC) 



Asinh(2Q;iAC)^ 



(27) 



where quantities ai are defined in Appendix A, see Eq.(A.12). Eq.(27) is one 
of the major results of the current paper: it generalises the results obtained 
by Mecozzi [10,11] and Turitsin et al. [25] for a single oscillator. As mentioned 
earlier, if the matrix coefficients diagonal, then each coefficient a, in 

(27) depends on Ui only, and propagator (27) factorises into the product of M 
identical functions, each having the form of Eq.(24). However, if the coupling 
matrix aij is nondiagonal (the case of nonlinearly coupled oscillators) , formula 
(27) cannot be reduced to a trivial product of identical components. 



4 Perturbational approach for zero-dispersion NLSE 



4-1 General remarks 



We are now in a position to build a general perturbation theory for the per- 
turbed zero-dispersion NLSE and its generalisations. In this section we will 
consider perturbations which do not contain a dispersive term. Any per- 
turbation operator N in Eq.(12) after regularisation yields a Fokker-Plank 
image A/ which will appear in the r.h.s. of Eq.(25), and will retain a small 
perturbation parameter e. Instead of equation (25) we arrive at the perturbed 
equation: 

^ + CkVk = sAfVK , (28) 

where K — M x N, M is the dimensionality of the base system and N is the 
number of samples. Since the coupling matrix a^j is block-diagonal, wc have 
Ck = J^iLi^^^M- When we consider local perturbations to system (12) (or 
regularised system (15)), the FP-image of the perturbation also decomposes 
into a sum of independent operators, J\f = J2iLiJ^i, each acting on a specific 
set of variables {r*, belonging to a separate block. Using propagator (27) 
we arrive at the following expansion (we write only the first corrections to the 
unperturbed propagator) : 

V^(q\ q'; = P^'iql 9°; C) + ^Af)(^| q'; + Ois') , (29) 
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where for P^^ we have: 



^(f)(g|9';C) = JdCjd-q'^gKiC-C; (t>K- 4>'K;rKyK) 

Q 

X N' \gK{C'\ <I>'k - <Pk; r^K, r'^)] , (30) 



AT 



Qk{C - C; <t>K - <t>'K\ r'^) = n - C; 0W - r-W, O , (3i) 

i=l 

oo 2iT oo 2lT 

J d~q'^ ^ J dr[ r[ J d<P[ ■ ■ ■ J dr'^ r'^ J d<P'^ . (32) 



K 

The prime in J\f' indicates that operator jQ" acts on the primed variables. 
Note that since perturbation enters only advection terms in Langevin equation 
(4) (or (15)) its Fokker-Planck image has a form of divergence: A/" = V ■ 
(...) (see equation (5)). Recalling also that the propagator Qm is normalised 
we can verify from equations (29)-(30) that the overall normalisation of the 
function Vk is preserved since jQdq'f^P^>^{q\ q'; () — 0. The PDF in cartesian 
coordinates V^{q\q°:, C) can be restored via the substitution 

V^'iqlq'; = V^'irix, y), ct>{x, y)\ r\x\ y'), ct>\x', /); C) 

In a similar fashion, continuing this recurrent procedure we can find the higher 
corrections to the PDF (29). 

Let us now consider explicitly some examples of nondispersive perturbations. 
4-2 High order nonlinearities 



First we consider the easiest case where the dimensionality of the correspond- 
ing base system is one and a perturbation is local. It seems quite natural 
to choose the higher nonlinear terms as a perturbation for the single zero- 
dispersion noisy NLSE (12): 

N[u]^ij2cku\u\^'+\ (33) 

k=l 

where q is any positive integer number and are some real phenomenological 
coefficients. Obviously, the condition Ni[u] = Ni[ui], is fulfilled: this pertur- 
bation is local. After the regularisation the FP-image of this operator appears 



15 



as a sum, jv = YnLi-^i^ with N being the number of samples and jVj given 
by: 

Up to the first order in e the solution of (28) is then given by Eq.(29), where 
one should substitute in Eq.(31) the Green function of the corresponding di- 
mensionality, Eq.(24): Qm — Gi- Any local perturbation can be handled in 
the same way. 

4.3 More general systems: two and more zero- dispersion NLSEs coupled via 
perturbations 

Now we proceed to more complicated systems, e.g. when we need to consider 
a set of equations (12) with coupling perturbations. In this subsection the 
number of dimensions of the corresponding base systems is still one, but the 
perturbation now cannot be treated as local. For instance if we consider two 
nonlinear couplers or a single birefringent fibre, then in place of Eq.(12) we 
obtain the following system: 

i|wp« + V(t,C)+eiViKC),t;(C)], 

(35) 

i\v\% + rf{t,0+eNMOAO]- 

Note that complex noises, rj^ and rj^^ are the independent WGNs, and may 
possibly have different intensities Di and Z^2- The perturbation operators, 
A^^[m,v] and (or) A^^[m,v], are supposed to intermix the components u and 
V because otherwise the latter would be uncorrelated and one could merely 
consider the first and the second equations of system (35) independently using 
the results of the previous subsection. 

It is convenient to introduce a second upper index to distinguish the vari- 
ables for u (superscript "1") and v (superscript "2"), and opt for the fol- 
lowing indexing rule: after the regularisation we introduce vectors r2N — 
(r},...,rf,r^, ...,r^), and </)2jv = {(t>\, ■ ■ ■ Ai Ah ■ ■ ■ ^2)^ where 

r\ = \u{i^t)l rl = \v{i^t)l i = l,...Ar, 
0i =Arg[«(iAt)], 4 = Arg[T;(iAt)], i = l,...7V. ^ > 

We seek the PDF of perturbed system, Eq.(28): VxivK, </>a-| r%, 0^; C), K = 
2N. The perturbation in (35) couples the dynamics (and hence the statistics) 
of variables {r\,(}){) and {r2,4>2) for i = 1, . . . , A^. Since we have only two 
(possibly) different noise intensities, Di and D2, for this case we can rewrite 



du 

dC 
dv 
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general expression (31) in more simple explicit form: 



N 
i=l 

xGi'iC-C; <P^-4>'i;4,r'^), (37) 



where the superscripts of the propagators mean that one ought to insert the 
corresponding noise intensity in the expression for Qi, Eq.(24). We can then 
employ the results of subsection 4.1. 

The generalisation for the case when wc have q zero-dispersion NLSEs coupled 
via the g-dimensional perturbation operator, 

is now straightforward. Firstly one should derive the FP-image Af for the per- 
turbation N. Next, by analogy with the case of two equations, one should 
introduce ii'-dimensional vectors {K = q x N): = [ri, ■ ■ ■ ,r/^) and (p^^ = 
{(pi, . . . , (px)- Again for the unperturbed system we have P^^Vk, (Pk\ q'k^ C) 
in the product form (20). To obtain the expansion for function we then use 
the formulae given in subsection 4.1. The expression for the propagator can 
be given in a simplified form similar to Eq.(37): it is a product of q equivalent 
multipliers Qi'' , k = l,...,q, where the corresponding noise intensities 
should be inserted in the expression for the propagator Eq.(24). 



4-3.1 Coupling perturbations of uncoupled systems 

Let us now provide some examples of coupling perturbations. 

• For birefringent fibres, with the coupling (~ s) between the different polar- 
isation components of either linearly or circularly polarised light, we have 
(see e.g. [14]): 



N[u, v] 



Ni[u,v 



Av\u\^ + Bu'^v* ^ 



(38) 



where A and B are the (real) constant coefficients of cross-phase modulation 

and four-waves mixing. The corresponding "polar" FP-image of in Eq. 
(28) by virtue of indexing scheme (36) can be presented as a sum A/" = 
Ef=iM, where A^ is, as usual, the number of time-samples, and for J\fi we 
have the following formula: 
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j^^ = -[A + B cos [2(0i - 4)]) 



J\2 



d 



d_ 
drl 



^' d<t>\ 

; d 



r+(ri)' 



9 



90: 



9rt 



(39) 



Now consider two nonlinear oscillators with small eigenfrequencies and weak 
linear coupling (both terms ~ e). This is a special case of self-trapping 
model [26] in the presence of the additive noise. This system is interesting 
in itself: it is integrable (in the noiseless case) and possesses a finite degree 
of freedom analogue of soliton solution (inhomogeneous state). Therefore 
this system can be reckoned as one of the simplest model examples, where 
we can investigate the action of noise on an integrable system. The coupling 
operator TV now is 



N[u, v] 



—UJU + ^{v — u) 



—CUV + 7(it ~ v) J 



(40) 



where uj and 7 are the real parameters. The number of independent variables 
is just four: u —>■ {ri, 4>i}, v {r2, 4>2}- The polar FP-image of N takes the 
form 



7COS(02 - 0i) I — TT— + 



ji d(pi r2 d(p2 J 
The expression for the propagator is very simple: 



(41) 



-Gi'iC-C; 0i-</''i;ri,ri)6?f^(C-C'; 02 - 0^;r2,r^) , 

where, as before, the expression for each Q^'' is given by Eq.(24). 
• Consider a circular nonlinear fibre array comprising q fibres with weak e) 
linear coupling. Then instead of (35) we acquire a set of q equations for q 
functions Vj: v — {vi, . . . ,Vq). The g'-dimensional coupling operator now 
becomes [14]: 



N[v] 



iC 



Vl + 



\ Vq^i +Vi J 



(42) 
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with C a real constant. Again it is convenient to introduce the double index- 
ation scheme: the upper index wiU indicate the sample number, the lower 
one corresponds to the field number: t'i(iAt) — >• {r\,(j)\}, . . . ,Vq{i^t) —>■ 
{r*,0*}. Then the image of the perturbation (42) can be presented as a 
sum, Jy — J2iLi-^i (with N the number of samples), where 



n=l 



Kirl<Pi) COS. 



_d_ 

dri 



sin (f)\ d 



( . j d cos( 



and we used denotations 



_d_ 



(43) 



KK. = r\-i sin + r;+i sin 0^^+^ , 
Bn{r% (jl) = -r;_i cos - r^+i cos 0: 



Here we have assumed that (rg, 0o) = (^g) 0g) and (r*_,_;^, 0^_,_;^) = (r},0*j^). 
The explicit form for the propagator is now: 



N 



Qk{C - C; 4>K - 4>'k; vk, r'^) = n 11 ^^(C - C; 0l - 0^; 4, ^10 , 

with K ^qx N. 



i=l k=l 



4-3.2 Birefringent fibre with weak four-wave mixing 

So far the original unperturbed system has always been equivalent to the 
system of uncoupled oscillators, i.e. matrix ajj in (8), (9) and (27) has always 
been an identity matrix and thus the dimensionality of the base system has 
been one. This is not the case for a birefringent fibre with a strong cross- 
phase modulation, i.e. if the inequality e holds in Eq.(38). Then one has 
to consider a more general perturbed system at zero dispersion: 

i{\u\' + A\v\')u + v\t, + eN'[u{C),v{C)] , 

(44) 

i{A\u\' + \v\')v + ri% + sN'HO, v(0] , 

where A is a real constant. Putting A = will correspond to system (35). We 
see that now the dimensionality of the base system is two, the matrix aij of 
the base system has the elements an = 022 — i, 0,12 — 021 = ^. As a simplest 
example of a perturbation for system (44) we choose the four-wave mixing (a 
special case of (38)): 

N[u,v]=iB\ , (45) 



du 

dC 
dv 
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with real constant B; the intensities of the noises, 1]^ and 77^, can, of course, 
be different: Di and UnUke the case of system (38), now the perturbation 
(45) is a local one: Ni[u, v] = Ni[ui, Vi]. The "polar" FP-image of N has been 
already derived as a part of expression (39), responsible for the four- wave 
mixing, i.e. Af = J2iLi M with: 



N/B^-cos 2(0^- 



„n2 



d 



J\2 



d 



+ r\ r\ sin 



X r 



d 



. d 



dr\ 

The propagator is given by 



2(0i - 0^2) 

(46) 



N 



Qk{C - C; <t>K - ^'k, rK, r'^) = n ^2(C - C; 4 - ; , (47) 



i=l 



with the corresponding noise intensities and coupling matrix elements in- 
serted into general definition (27). 

The cases of coupling (nonlocal) perturbations and the higher number of di- 
mensions of the corresponding base system can be treated analogously. 



4-4 Discrete self-trapping equation 



It is instructive to give explicit expressions for noisy DSTE. In this system 
the discreteness is postulated by the model itself, the dimensionality of the 
base system is one, and we can readily find the perturbation of FP operator 
applying the results above. Now the perturbation reads (see equation (19)): 



N 



(48) 



where N defines the number of discrete sites in the chain (lattice). In this case 
the nonlocal perturbation (48) might couple all N oscillators. Its FP-image is 



AT 



n,j=l 



sm 



((f) j - (f)n) 



_d_ 

dr„ 



cos 



d 



(49) 



and the propagator of the system factorises into a product of N functions 
Qi% with the corresponding noise intensities. We can then apply the general 
formulae from subsection 4.1 to find the perturbed PDF. 
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5 Second order dispersion as a perturbation 



In this section we will consider more complicated examples of nonlocal per- 
turbations. We will concentrate on one particular type of perturbations which 

is nevertheless of great practical interest: the dispersive terms in Kerr sys- 
tems. Such terms can account for, for instance, a small residual dispersion 
in a dispersion-managed fibre link. To start we first study a scalar case of a 
single NLSE, and then proceed to vector generalisations, i.e. to the system of 
coupled NLSEs. 



5. 1 Scalar NLSE with weak second order dispersion. 



In this subsection we will consider the perturbation in terms of weak second 
order dispersion (SOD) in (12), with the dimensionality of the base system 
being one: 

(50) 

Here d = ±1 specifies the type of the dispersion (anomalous or normal). 
For the sake of simplicity we consider the case d = 1 since the choice of 
d = —1 only alters the sign in the appropriate formulae. As we work with 
discretely sampled data we need to write down a discrete analogue of the 
second derivative Eq.(50). To do so we first define the forward finite difference 
as Sfn = fn+i — fn- Thc arbitrary precision expansion for the ID second 
derivative, d^u{t)/dt^, reads as [27]: 



\n\l + d)u^ = -— J2 bn S^'ui , (51) 



where we have applied regularisation procedure (14) and 6" is the forward 
difference of the nth order: 



d'^Ui ^ Ui+i - 2ui + Ui-i, d^Ui ^ Ui+2 - ^Ui+i + 3ui - Ui-i , etc. 

The squared logarithm operator in representation (51) should be treated as 
a Taylor expansion in powers of S, and the expansion coefficients 6„ can be 
easily calculated: 62 = 1, ^'3 = —1, ^4 — 11/12, etc. 

Taking into account Eq.(51) we now write the discretised perturbation iVj as: 

00 

{AtfNM = iln2(l + d)u, = i Yl bnS^'u, (52) 

n=2 

Further on it is convenient to use the explicit expression for the difference 
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operators in terms of binomial decomposition [27]: 



<5"«i = E(-i 



j=0 




i+n-j 



The Fokker-Planck image of the perturbation J\f then becomes: 



Af oo n / \ 

(At)W=EEE(-i)C^^ 

i=ln=2j=0 \J / 

COS {^(f), 



^i+n-j 



n I i+n—j 



sm 



i+n-j 



•)- 



d 
d(t)i 



(53) 



Note that we write an infinite upper hmit of summation over the index n. 
However, of course, each difference is hmited by the total number of samples 
N , and all the variables outside the sampling interval should be discarded. 
The form of the perturbative expansion for the full PDF does not differ from 
that given in previous sections for the non-dispersive perturbations. The total 
propagator has the form of the A^-product of functions Qi, each having the 
same noise intensity D. 

The same procedure may be easily developed if one wishes to account for a 
perturbative term in the form of a derivative of an arbitrary order with either 
real or complex coefficient (~ e). For example, the dispersion of arbitrary 
order i^d^u/dfJ' would result in the change of the power of the logarithm 
operator in expressions (51), (52). Eventually one should arrive at the same 
formulae with a different set of coefficients which, in that case, have to 
be the Taylor coefficients for the expansion of the corresponding power of the 
logarithm operator (one will have to change the lower summation limit for n 
from 2 to A; in Eqs. (51-53) as well). It is worth mentioning that system (15) 
with the discrete term, given by (52) but truncated at n = 2, corresponds to 
the weakly discrete NLSE [3] with additive noise. The increase of the number 
of terms in expansion (52) (and consequently in (53)) obviously increases the 
precision at the expense of the computational time. 



In the same manner we could redefine the dispersion operator via the backward 
differences (see [27]) and obtain an expression very similar to (53). If we wish 
to preserve the time-inversion symmetry, we might use the symmetrised form 

instead of Eq.(51): A/" = {Sfforward + -^backward) f'^- 
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5.2 Manakov equations with weak second order dispersion 



Next we develop a perturbation theory accounting for the weak dispersion 
in vector systems, such as Manakov equations (17) (i.e. for the regularised 
system (18)). Now the number of dimensions of the base system is two and 
we consider the nonlocal perturbation of the form: 



We do not present here the explicit expressions for the PDF expansion as 
they are very similar to those given in the previous subsection, but outline 
only the general procedure. First we choose the same regularisation we used in 
subsection 4.3 for two uncoupled equations. Pull coupling matrix Uij comprise 
a set of 2 X 2 blocks placed on the main diagonal, with all elements within each 
block equal to 1. The expression for the Fokker-Plank image of perturbation 
(54) is very similar to Eq.(53). The only difference is that now we again have 
to endow the variables with an additional index to distinguish between the 
fields: Ui {r\,(p\), Vi — > {rl^cf)^). Now the expression for the FP-image of 
perturbation (54) consists of the two summands, J\fi and A2, each having 
the form of Eq.(53): the variables in Eq.(53) for A/i should bear additional 
index "1", and the same applies to A/2. The propagator for the entire sampled 
system is a direct product of with the corresponding (possibly different) 
noise intensities, Di and D2, and the elements of the coupling matrix for the 
base system, an — ai2 — 021 = 022 = 1- All statements from the previous 
subsection, concerning accounting for either the time inversion symmetry or 
a higher order dispersion apply to this case as well. 

Remark on higher systems of coupled NLSEs. Following the preceding exam- 
ple one can build perturbation expressions for the systems comprising more 
than two nonlinearly coupled NLSEs. The number of coupled NLSEs, say M, 
simply defines the dimensionality of the base oscillator system. If we have M 
coupled fields, Vi, . . . , vm, then introducing an additional index labelling the 
field number one finally arrives at the expression for the FP-image of a dis- 
persive perturbation in the form of a sum, J\f = J2k-^k, where each A/fc has 
the form of Eq.(53). The propagator for the whole sampled system is again a 
direct A^-product of Q^'' . 



6 Conclusion 

In the current paper we have proposed an approach for studying the stochastic 
dynamics of a noise-driven systems with Kerr- type nonlinearity. We presented 




(54) 
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the exact expression for the PDF (and the propagator) for a system of non- 
hnearly coupled oscillators. Aside from the study of nonlinear lattice models 
with the additive WGN, we showed how to apply these results to the descrip- 
tion of signal propagation in a weakly nonlinear dispersive media driven by 
WGN. We used the fact that in a system at zero dispersion the PDF for the 
output signal (field) can be obtained analytically. Using the propagators for a 
zero dispersion system we have built a perturbation theory and obtained the 
corrections to the unperturbed PDF in a variety of physically meaningful sit- 
uations. In particular we were able to obtain the corrections to the PDF in the 
presence of the nonlocal perturbation, such as a second order dispersion. The 
knowledge of the output signal statistics allows one to calculate a multitude 
of important quantities like the probability of error in the fibre optics commu- 
nication, or, in principle, Shannon capacity of the nonlinear communication 
channel, modelled by Eq.(12). 

Using the discrete time picture one can also examine the statistics of more 
complex systems than those, described in the current paper, like, for instance, 
a system of weakly coupled NLSE with (also weak) linear or nonlinear nonlocal 
perturbations (e.g. nonlinear dispersion). This would simply require increasing 
the number of dimensions in the corresponding PDF. It is also possible to 
consider the higher- dimensional noisy NLSE in the same fashion, e.g. the 2-1-1 
dimensional NLSE where the (discretised) weak second order derivatives are 
reckoned as a perturbation. 

We recognise that for the case of nonlocal perturbations, like higher order dis- 
persion, when the perturbation operator in the FPE becomes quite involved, 
the corrections to the PDF can only be calculated numerically. Let us recall 
that the traditional way of calculating the oiitput signal statistics numeri- 
cally is to use direct Monte Carlo simulations. The latter method, however 
has a number of significant drawbacks. Normally the calculations are quite 
time consuming and cannot be used for estimation of the tails of the PDFs, 
something which is crucial for the performance assessment in the fibre optical 
communications. Therefore the perturbation technique developed here, with 
the iteration procedure similar to (30) performed numerically, may present an 
important method of choice for the systems where the direct Monte Carlo ap- 
proach fails. The proposed method relies only on the approximate evaluation 
of integrals and infinite sums without direct integration of the noisy equation. 
And finally, the pertrubative approach devised in the present paper does not 
require the noise to be small - which is usually implied in the majority of 
the approximate methods in stochastic dynamics. Therefore it can be used to 
analyse systems far from the equilibrium where the impact of noise on the 
signal degrees of freedom becomes significant. 

In conclusion we also point out that the exact analytical results concerning 
the statistics of nonhnear coupled systems are still few and far between. Thus 
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we think that the exact general expression for the propagator of coupled noisy 
nonlinear oscillators (27) (which coincides with a retarded Green function), is 
quite important in itself and can serve as a "reference point" for the studies 
of similar systems. 
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A Eigenfunctions of the Fokker-Planck operator 

In this section we calculate the eigenfunctions and eigenvalues of the FPE for 
single oscillator (22) and for system of nonlinearly coupled oscillators (25). We 
start from the former case. We wish to solve the eigenvalue problem: 



Since the obvious periodical boundary conditions in it is convenient to per- 
form a Fourier transform with respect to 0: 



(A.l) 



oo 



^(r, A0) 



J2 e''^^^*.(r) 




(A.2) 







After the transform the operator Ci reads as 




(A.3) 



To solve the eigenfunction equation we make the following substitution: 



*,(y) = y 1^^1/2 exp(-a,y/2)(^(y) , (A.4) 
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with — {s — 2 1 + \u\ a„)/4:. Equation (A. 5) has the canonical form of the 
degenerated hypergeometric equation [28]. Since we want our eigenfunctions 
to decrease at infinity we must demand that Zgi^/au = n, n = 0,1. . . .. This 
will determine the discrete eigenvalues Sm, of (A.l), while the corresponding 
eigenfunctions are expressed via generalised Laguerre functions L^(. . .) 
as [28]: 

Sni, = 2Q;,.(2n + 1 + |i/|) , 

1/2 / f \ ^/"^ 

*„.(r,0)=e-^(^) ((^^^j .H/2exp[-./2]LH(,), (A.6) 
z — ai,r^ /D . 

Operator Ci is of course not Hermitian. Therefore we need to introduce a set 
of left-eigenfunctions "^nvi which are the eigenfunctions of the adjoint operator 
C[ with eigenvalues s„i,: 



«»„(r-,0)=e'-*(^)'" (^-^y'%.M/^exp[-.721L!rl(..). 
Eigenfunctions (A. 6), (A. 7) form a biorthogonal system, i.e. 

27r oo 

j d(t) j drr^*^,{r, (l))^n'.'{r, 0) = S^n' ^v.' ■ (A.8) 



Making use of the identity (see [28]): 
> n! ^ ^ = L exp(-w ) 1^ 2- , 



|w| < 1 , 



(A.9) 

analytically extended to the complex plane, in the limit to ^ 1 we can obtain 
the closure relation for the left and right eigenfunctions: 

^ ^;,(r, 0)^„.(r', 00 = 5{r - r') 5{4> - 4>') (A.IO) 

n,v 

(factor in the r.h.s. is merely a jacobian). 

In the same way, generalising the derivation of the elementary basis (A.6), 
(A. 7), we can built left- and right-eigenfunction for the general 2M-dimensional 
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operator Lm in (25). Using the Dirac notations one gains 



M 



ni, . . . ,nM,Jyi, ■ ■ ■ ,J^m) = 



TT 



M/2 



1=1 



X 



rii, . . . , um, t^i, ■ ■ ■ , t^m\ 



TT 



M/2 



M 



i=l 





riil 


/A 


_{rii + 



1/2 



1/2 



x<i-i/V</^Li^i(.;), 

Here the quantities Q;i(z/i, . . . , z/m) are defined as: 

1/2 



(1 + i) (f Ef=i a, for E^li > , 

. (1 - i) (f |Elii a,,z/,|) for Efii < . 



The corresponding (right) eigenvalues of operator Cm are 



M 



'ni,...,nMfi,---,i'M 



2^a,(2ni + l + |z/i|). 
1=1 



(A.ll) 



(A.12) 



(A.13) 



Since wc know the exphcit expressions for eigenvahics, we know the damping 
eigenratcs (7 = |Re[s]|) for each eigenmode of FPE. These may be important 
for the investigation of the long scale (large values of Q evolution of the 
system statistics, because for the case ( ^ 00 the statistics of the system is 
determined by a single eigenmode with the lowest value of 7. We can see that 
the damping rate for each eigenmode is proportional to 7(2/) ~ D^/^i/^/^ for a 
single oscillator or, for a particular realisation (z/i, . . . , z/^r), to 7(2/1, . . . , Pn) ~ 
Yl!ij=iDl^'^Vi{ajiVjY/'^ for a system of nonlinearly coupled oscillators. Note 
that the real parts of eigenvalues (A. 6), (A. 13) are always positive which means 
that equations (22), (25) do not possess a stationary solution. 
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